

KS) 

( *** C -T^) 


Model Calculations for Diffuse Molecular Clouds 

by 

A. E. Glassgold and William D. Longer 

New York University 
New York, New York 10003 


Received 


** Uaiy i 

... CSc £ 20 H 


# 74 - 


22 3 u 5 




- 2 - 


PRRCRn,i>j Cr 


P.A a m 


m ' ANK NOT FIT.MRD 

ABSTRACT 


A steady state Isobar ic cloud model is developed. The -pressure, thermal, 
electrical, and chemical .balance equations are solved simultaneously with a 

simple one -dimensional approximation to the equation of radiative transfer 
appropriate to diffuse clouds. Cooling is mainly by CII fine structure 
transitions, and a variety of heating mechanisms are considered. Particular 
attention is given to the abundance variation of fL using an extension of 

ft 

the work of Hollenhach, V/erner, an.d Saipeter. Inhomogeneous density distri- 
butions are obtained because of the attenuation of the interstellar l JV 
field and the conversion from atomic to molecular hydrogen . The effects 
of changing the mode], parameters are described and the applicability of the 
model to OAO-3 observations is discussed. Good qualitative agreement with 
the fractional K 0 abundance determinations has been obtained. The observed 
kinetic temperatures near SQ C K can also be achieved by grain photoelectron 
heating. The problem of the electron density is solved taking special account 
of the various hydrogen ions as veil as heavier ones. 

CII column density determinations are also emphasized. 


The roles of Cl and 
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I . Introduce i on_ 

Observational information on interstellar clouds is rapidly increasing 
to such an extent that significant tests of theoretical models now seen feasible. 
We have recently presented a preliminary discussion (Glassgold and Longer , 1973c) 
of the H and G abundances measured for 15 or so diffuse clouds with the Cop er- 
nicus 0A0-3 spacecraft (Spitzer et al, 1973; Morton et al, 1973). This article 
is the first relatively complete report of our work, and Includes a full dis- 
cussion of cirri properties and parameters as they relate to the 0A0-3 results. 

In future publications we will discuss thicker clouds and other molecules 
besides 

Our present model is basically a local steady sir '-e model in which all 

quantities vary iron point to point except for the pressure (p^), which is 

uniform. The inhomogeneities arise from the attenuation of the external 

radiations which heat, ionize, and dissociate the cloud constituents. The 

most import arn external radiation is the Iff component of the interstellar 

radiation field. No doubt a steady state model is over idealized, but at 

least it permits a beginning to understanding how the many physical phenomena 

interact in interstellar cloud material. The isobaric condition should be 

appropriate if the sound crossing time is less than the most lifetime of the 

6 7 

cloud. For the QAO-3 clouds, the crossing tine lies in the range 10 - 10 

years . 


The emphasis in this first report is on those components of an inter- 
stellar cloud which are most important for determining its thermal properties, 
i.e. the gas temperature T and density in Thus we include the balance between 
R and because hydrogen is the most abundant species, and the balance 
between Cl and CII because carbon is the most important cooling agent at 
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the temperatures characteristic of diffuse clouds l®3 K). The inter- 
stellar grains are also included because they significantly attenuate the 
interstellar UV radiation field, catalyze the formation! of and possibly 
heat the gas via photoelectron emission. The electron (density n is also 
calculated because the electrons play a basic role in imany of the underlying 
atomic reactions. The chemical equilibrium for caleiuim is included since 

it is relevant to determinations of n from measurements of the column =r; 

e 

densities of Cal and Call (White, 1973 and references elted therein). 

Other constituents are not discussed here in detail, although our program 
has been constructed to deal with them. We have not attempted to analyze 
and to fit exhaustively the observations for any one particular cloud, 
but instead have concentrated on the general properties of the average 
cloud, represented by the ensemble of 0A0-3 cloud observations. 1*. 

One of our most important concerns is the heating of the clouds,: . ■ 

■ since the cooling mechanisms for diffuse clouds seem tea be fairly well 
established (Dalgamo and McCray, 1972). The problem Iteere is to achieve 
relatively high temperatures in the neighborhood of IOgP K (Hughes et al, 
1971; Radhakrishnan et al, 1972; and Spitzer and Cochran, 1973). Although 
the observed C depletion (Merton et al, 1973) helps, significant grain 
photoelectron heating seems to be required, as advocated by Watson (1972). 

Earlier work most closely related to this investigation is that of 
Hollenbach, Werner, and Salpeter (1971, referred to henceforth as HWS), who 
demonstrated quantitatively that will be concentrated in the interior of 
clouds by self- shielding of the dissociating UV radiation. HWS 
considered only the balance between UV dissociation of through the Lyman 
band lines (Solomon and Wickramasinglie, 1969) and formation on grains 
(Hollenbach and Salpeter, 1970) for constant total hydrogen density n and 
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temperature T. In this work we include the significant thermodynamic 
consequences of this balance and of the associated attenuation of the 
interstellar radiation field, as well as other chemical balances involving 
carbon* electrons, hydrogen ions, etc. The attenuation leads to significant 
inhomogeneities especially in n, which then feed back to affect the con- 
centration. 

In Section II we outline our calculations, and in Section III we >' 
describe the results obtained, which include the effect of varying the most 
important parameters and suggested ranges for them. A Tbrief summary and 
some conclusions are given in Section IV. 

II. Calculation 

A. Overall Structure of the Calculation 
We adopt a one- dimensional steady state model isa which the 
density and column density of species X at position x ame denoted 
respectively by n(X,x) and 

N(X,x) = C dx T n(X,x* ), (!) 

0 

; 

x - 0 being the cloud surface. However, we usually suppress the 
x-dependence with the understanding that everything except pressure 
depends on x. In place of x we often use the total column density 
N of hydrogen to measure distance into the cloud 

N = N(H) + 2N(H 2 ) (2) 

The radiation field intensity is I(A,x). Although the calculation 
involves a semi -infinite cloud medium, the treatment of the radiation 
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transfer discussed below is done so that finite-cloud effects are 
included in an approximate way. 

The basic equations in the model are steady-state balances. 


pressure: 

^n^(x) + 2 n(X,x) J kT s p^ 
X 

(3) 

thermal: 

A(x) * r(x) 

U) 

electrical: 

n (x) “ 2 n (X,x) 

x i 

(5) 

chemical: 

P(X,x) * D(X,x) 

(6) 

In these equations, 

P(x) and D(x) are respectively the production 


and destruction rates for X, n^(X,x) is the total ionic density for 

species X, and A and V are respectively the total cooling and 

heating rates. Effects of macroscopic motion and gradients are 

neglected. The local nature of these steady state balance equations 

has been emphasized by Eqs. ( 3 )— ( 6 ) . Each quantity also depends in 

principle on temperature T, on densities n(X), and on the radiation 

field I. Once I is known, there are just enough equations to solve 

for . all densities n(X) and the temperature as a function of x and 

the pressure p ♦ (If there are s species, there are s + 2 unknowns 

including n and T, and Eqs, (3)-( 6) provide s + 2 equations since 
0 

(6) is equivalent to s homogeneous equations.) Of course the radiation 
field If Vx) depends on the solutions to ( 3 )-( 6 ) for all points x* 
such that 0 $. x 1 <. x. 


Thus the solutions for the density and 
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temperature profiles n(X,x) n-C'x), T(x) and the intensity I(*,x) must 
be obtained in a self-consistent manner. 

B. Restricted Set of Balance Equations for Diffuse Clouds 

By diffuse clouds we mean those with column densities to the 

21-2 < 
center N < 10 cm or with visual extinction 'u 1. The 0A0-3 

clouds fulfill this condition. Cosmic ray protons with energies as 

small as 2 MeV will hardly be attenuated by such clouds. On the 

other hand, the incident radiation field is significantly attenuated 

by the cloud material, especially in the far UV by the dust grains lhr 

and molecules. As long as the temperature T 5 100°K, CII is the 

most effective of the atomic cooling agents. (The abundance of Cl is 

too small to contribute significantly to cooling. ) As for molecules, 

is not very effective (Glassgold and Larger, 1973a )> and neither 

CO or HD are observed in sufficient abundance (Jenkins et al, 1973$ 

Spitzer et al, 1973) to be competitive with CII (Glassgold and Langer,*-tb 

1973a; Dalgarno and Wright, 1973). Thus the atomic species which 

determines the cloud structure are H, Cl, CII, grains, and electrons. 

We now discuss their balance equations and those for Ca, which serves rirl 4 

as a useful prototype for further applications of our model. 

We do not consider the gain or loss of grain atoms, but simply 
specify the average grain area per gas atom inferred from the 0A0-2 
and 3 extinction measurements (Bless and Savage, 1973; York et al, 1973 )* 
l * e • , 

Sgr ^ °gr ~ 1.036 N 


(7) 
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Here 0 (X) is the mean total cross section of a [grain for radiation 

gr 

of wavelength X , and E (X) = n (X)/n is the corresponding grain 

gr gr 

to gas density ratio, ignoring He and heavy atoms. In interpreting 

the extinction measurements, we replace the column density ratio 

N /N by K . Eq. (7) enters into the attenuation of the incident 
gr gr 

radiation and the catalysis of formation, both (discussed in more - - 
detail below* , * 

The formation and destruction of is treated along the lines of 
HWS, but cosmic ray destruction (Solomon and Werner, 1971) and the effects 
of ion molecule reactions are also included: 


. ■ t \ ( <ato(H_) i 

R^OOn = ^ 2 (H 2 ) ♦ n(H 2 ,I) ) n(H 2 ) - ■) 


( 8 ) 


ion react 


On the left .side we have production on grains with, sa rate constant 

n 


' R 1 = 1/2 Y v th (n) -E- A 


gr 


(9) 


where y is an efficiency factor ( ^ 0.3 for T = 1Q1D° K, Hollenbach ' 

and Salpeter, 1970), v^fll) is " blie mean s P eed of H atoms, 

n is the total density of grains and their Mean area for catalyzing 

gr gr 

formation. By adopting a ,, typical ,, grain model with 

op o -17 1 -1 

n A /n * 3,8 x 10 cm , HWS estimated R_ * 10 x cm^s , fairly 
gr gr l 

independent of temperature for T < 100° K. Much earlier, McCrea 

—17 3 -1 

and McNally (i960) had estimated - 7 x 10 cmr s by using 
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somewhat smaller and denser grains and the maximum efficiency y * 1. 

On the other hand, the large far UV-extinction observed by 0A0-2 and -3 

imply from (7) values for E a as large as 3 x 10 cm for the 

gr gr 

. -21 2 
typical reddened star and 5 x 10 cm for £ Oph. These measurements 

seem to suggest a major component in the grain distribution (York et 

al, 1973) of small particles of linear dimension ^ 150 A . If the 

corresponding area is assumed to be about as effective as the larger 

grains considered by HWS; then their estimate of R^ should be revised 

upward to 'v 5 x 10 ^cm^s ^ for the typical reddened star (Spitzer 

-17 3 —1 

and Cochran, 1973) and to 7 x 10 cms for £ Oph. Because of the 

considerable uncertainty in these estimates, we have investigated a 
range for R^ from 5 x 10 - lO'^^cm^s*"^. 

The right side of (8) contains the ionization rate by cosmic 
rays and the photo- dissociation rate p. Although direct dissociation 
of by protons with cosmic ray energies is unimportant, dissociation 
by primary electrons is not (Glassgold and Langer, 1973b). Close to 
98% of the inelastic events produce + e and 2% H + H + e (Adamcyzk 

et al, 1966). These reactions are summarizrd in the table of Appendix 
A. Because the low energy cosmic ray flux is not directly measured, 
we have considered a range for r (H), the primary ionization rate 

}r 

of atomic hydrogen, from 10~ s - 10 ? s . Photodissociation 

of below 13.6 eV requires a two-step process in which discrete 

line excitation is followed by dissociation (Solomon in Field et al 

1966). The rate r^H^jl) will depend on the values of I (X,x) at 

00 

appropriate wave lengths in the range from 1108 A - 912 A , which 
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depend sensitively on position x. Our calculation is based on the 

method given in the appendix of HWS, and will be discussed briefly 

later. For a very thin cloud and Habing's (1963) field, we obtain 
-10 -1 

t) = 0.44 x 10 s . 
o 

Equation (3) also contains a term arising from ion-molecule, 
charge-exchange, and dissociative recombination reactions of the 

■ t “V* 

cosmic ray produced ions and . These reactions are important 
for heating and ionizing molecular clouds (Glassgold and Langer, 
1973a, b), and for the production of heavier molecules (Watson, 1973; 
Herbst and Klemperer, 1973; Dalgamo, Oppenheimer, and Berry, 1973). 


The term 




ion react 


in (8) is of the order of magnitude 


of and becomes important only when n £ 5 . These reactions 
are incorporated in our program, and we discuss them below in Appendix. 1 
' A in connection with the electrical balance equations. The gas density 
n, excluding He and heavy atoms, is the hydrogen mass density in units 

of v 


n = n(H) + 2n(H 2 ). 


( 10 ) 


When 


/ dn(H )\ 

I -7T — I and £ 0 n(H 0 ) are ignored, the fraction 

x / ion react d * 


of hydrogen mass in molecular form is 


2n(H_) 


f = 


n 


= (1 + 


n 


-1 


R i n 


(id 


The chemical balances for C, Ca and other heavy atoms are 
determined by the familiar relations (Spitzer, 1968); 
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n(X r + x ) G(X r ) ( 12 ) 

“< X r + iK 

where a is the radiative recombination rate coefficient and G is the 
photo ionization rate 

VV 

G(X r ,x) = $ o dX cr.'X r ,A) j c 0Qnt (xV K*,*) J - ( 13 ) 

912 A L 

Cosmic ray ionization has been ignored. In ( 13 ) , a ^ is the photo 
ionization cross section and c is a renormalization factor converting 
intensity to flux discussed below. For Cl and Call, may be replaced 
by their threshold values, whereas for Cal the complete observed cross 
section (Carter et al, 1971) increased by an overall factor of 2 
(Mcllrath and Sandeman, 1972) is used. The values needed to evaluate 
(12) for n( CII )/n(CI ), n( Call )/n( Cal ), and n( Calll )/n( Call ) are given 
in Table 1, as are values of G/a at zero optical depth for a very 
diffuse cloud. The reduction in G/a by UV attenuation is included 
in our calculations; the method for treatin' 3 ’ the attenuation is 
discussed below. In addition to (12), we have the total abundance 
of species X 

n(x) = Z n(X j , (14) 

r 

and the relative abundance 


£(X) = n(X)/n 


(15) 
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We consider only 01, Oil, Cal, Call, Calll because A^(CIII), A^(Cal7) < 912® . 

In accord with the preliminary 0A0-3 C-depletion determinations (Morton 

-5 _/ 

et al, 1973)? we consider values for £ (c) in the range from 10 - 10 . 

Although Ca is generally under abundant, we do not concern ourselves 

with choices of £ (Ca) since Calll (and thus the total Ca abundance) 

is not measured; also Ca does not affect the thermodynamic structure 

of the cloud. «■ 

The ion densities required to evaluate the electron density 

from (5) may be classified into two main groups, those of the light 

atoms (H and He) and those of the heavy atoms (0, C, N, etc.). In 

our model, the ionizing agents are low energy cosmic rays and the ~ r ‘"' 'V* 

stellar radiation field cut-off short of 912 £, This type of 

radiation field can not ionize the light atoms, but cosmic rays 

^ will • with characteristic rates which we take to be in the range 
—18 —1 — 16 — D 

10 s - 10 s The radiation field ionizes heavy atoms with e.i.u .’.or:."? 

ionization potentials I < 13.59& eV with a characteristic rate' 

-10 -1 

F i * 10 s for diffuse clouds. Despite significant attenuation 
in the far UV, photo ionization is the l :>st important ionization .-»<■’ . n .. 
mechanism in our model for such heavy atoms as C, Mg, Si, etc., in - 
gross qualitative agreement with 0A0-3 cloud abundance determinations 
(Morton et al, 1973). * , 

All the ions recombine radiatively with characteristic rate 

coefficient a - 7 x 10 ^ (T/100) cm^s ^ except for the 

t ( + 

molecular ions. The cosmic ray produced ions are rapidly 

•f 

transformed into ions and then recombine dissociatively, i.e.. 
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p + H^ + + e + p f 


(l6.a) 

H 2 + + H 2 -*■ H + + H 


(16. b) 


, ' \ 


+ 



+ e + H 2 + H and 3H 


(l6.c ) 


The detailed analysis of ion densities in Appendix a shows n(H 2 + ), 

) « n^ . This lo to be expected from the very large rate 

coefficients for (I6.b): K 2 = 2 x lO^cm^s"" 1 (Trujillo and Neynaber, 

1968; Bowers et al, 1969) and for (l6.c): ^ = 4 x 10~^ T^cm^ 1 

(Leu et al, 1973). The balance equations for the heavy atoms were 

already given in (12), and a complete study of the ionization balance 

is carried out in Appendix A. In our program we use the results( A-30 )- 

(A-40), which apply if N < 10 cm" 4 and x ^ 0.1. Our result for the 

e 

electron concentration becomes relatively simple if n(CIl) » n(Cl): 

X e = V2 (Sj^ + /s 2 ^ + 4 C'/n a (H*) ), (17) 

where ^HA is the abundance of ionized heavy elements. 


P HA 


Z g(X) 
X 

I(X) < 1(H) 


and 
C’ = 


j^l + 1,/2 ? 2^ H 2 ^ 


V 


+ (K 2 /2 - K^f 


(l-f) 


J. 


* 1/2 c(H + )f + £(He) c(He) a ^ H j 

a(He ) 


( 18 ) 


(19) 


is an effective cosmic ray ionization rate. The reaction rates K and K 
are given in Table A and £(He) » 1/14 (Cameron, 1973); also a(He + ) ~ a(R 
(Seaton, 1951). 


CM + 
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The dominant contribution to (18) is expected to come from C. Most 
heavy elements are observed to be significantly depleted in the 0A0-3 
clouds, but there is some uncertainty whether £(N) should be included 


in (IS) (Morton et al, 1973). The maximum value of obtained 
from solar abundances is 6 x 10 but the observed depletion (Morton 
et al, 1973) suggests the range 5 x 10 ^ - 1 x 10~^, with the smaller 


number representing overall depletion by 10 and no Nil. For 0A0-3 

-12 3 -1 

clouds, we might use the parameters (1 - f = 1, a * 6.4 x 10 cm s , 

_Q 

n » 25 cm ) to characterize the surface region, and (l - f B 0.2, 

-12 3-1 -3 

a, ® 8,3 x 10 cm s ,n = 50 cm )for the interior of the largest clouds;. 

the corresponding values of x are respectively 8 x 10^ and 

2.5 x 10 ** when £ (H) * 10 . As a result or depletion, cosmic 

P 

-''rays will be the dominant ionization mechanism in 'diffuse clouds even 

-17 -1 

if r (H)is as small as 10 s 

. p 

C. The Interstellar UV Field and Its Attenuation 

( U M 1 I I M ' — » » 11 I 

The thermodynamic structure of interstellar clouds is strongly 

o 

affected by the interstellar radiation field in the band from 912 A 
o 

to 1101 A, the ionization threshold for Cl. The I^man band lines 
o 

begin at 1108 A, although the ones important for dissociation do not 
o 

occur until 1026 A (\> f = 6 v = 0). The radiation field down to 


2028 £ is also required for the chemical balance of Ca. 

In discussing the intensity I(X ,x) as a function of wavelength 


X, we ignore reradiation by both grains and molecules at wavelengths 
X * > X . This would be a serious error if we were treating the 
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population distribution of the molecules or the thermal properties cf 
the grains, but we believe that such effects are not important in 
considering the total abundance . This is in accord with earlier 
work by MS; Black and Dalgamo (1973) have recently considered the 
effects of the fluorescent radiation on the population of the 
molecules. Elastic scattering by molecules is ignored, whereas 
elastic scattering by grains is assumed to occur mainly in the near 
forward directions. We thus adopt a transfer model in which each 
wavelength A is individually attenuated by the Iyman band lines 
and the grains. Because the line absorption is fundamentally different 
from the grain continuum absorption, we have two distinct phenomena to 
discuss. 

There are other sources of line and continuum absorption, but- 

the 60 Lyman lines (v 1 =0-19 **- v = 0) and the grains are the 

most important in the spectral range of interest. Because of the low. 

abundance of heavy atoms and the much higher- frequency absorption 

lines of He, the only real contenders for overlapping lines are the 

Lyman series of H and the Werner bands of The probability of :> ’ 

overlap is increased by the fact that considerable absorption is 

occurring in the radiative wings of the lines, i.e. the lines 

have an effective width AvV{fN ) ? where f is the oscillator strength. 

o 

We find that L 0 and L shield the R y 7i) - 1024.96 A -and 
p Y 6 +■ 0 

o 

H^i ^ q( 0) " 971.99 A transitions, respectively. (We use the 

standard notation for the X ->■ B and B X Lyman band transitions: 
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v'J* +■ vJ and v'J' + v' 1 J* 1 or K M , depending on whether dissociation 

occurs or not. ) Momentarily assuming a constant UV flux in the 1025 
o 

925 A band where most of the dissociation occurs, we estimate that 
these two lines would contribute only 2.7% of the dissociations. 

Because of the reduced oscillator strengths for both the Lyman series 
and Lyman band lines near the former's limit, shielding of the upward 
transitions \>‘ *■ o for v’ £ 19 is probably not effective for 

0A0-3 clouds. Moreover, the three components of 19 0 , which are 

nearly coincident with Lyman lines near the series limit, contribute 
only 1% of the (dissociations. The predicted decrease in UV flux in 
this region (Jura, 1974) also implies a further reduction in the 
role of these lines. Thus the overlapping Lyman series reduces the 
dissociation rate by only ^ 3%, which is negligible. Overlap with 

S' 

the Werner bands also appears to be unimportant in this context, 

there being only one hear coincidence. The Werner bands 

o 

accessible with radiation from 912 - 1009 A have negligible 

dissociation probabilities (Stephens and Dalgarno, 1972). The 

photoionization of heavy atoms does not appear to be an important 

o 

attenuator of the UV continuum below 2000 A because of their depletion 

and the large grain cross section associated with the 1/X extinction. 

The most abundant heavy atoms with I < 13.6 eV, Mg, Si, Fe, and S, 

have a total solar abundance of ^ 10 ^ and a very uncertain average 

-17 2 

photoionization cross section ^10 cm . Assuming depletion 

by 'v 10 \ they contribute an attenuation cross section of 

-22 2 -21 2 
^ 10 cm per gas atom as compared with 2 x 10 cm from grains. 


>o 
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As mentioned earlier, our treatment of the Lyrnaai band line 
photo dissociation of follows closely that described in the 
appendix of HWS, so we only mention a few important changes. 

As in HWS, we assume that the molecules are in the ground 
electronic state with v « 0, J = 0 and 1. This low state of 
excitation seems to apply to most but not all of the: (QA0-3 
observations (Spitzer and Cochran, 1973). We have included the 
possibility of arbitrary populations Pj for the two ground states, 
but have found the calculated abundance of H 2 to be insensitive to 
whether Pj is given by the high temperature limit P_,/)P « 3 used 

by HWS or by the observations of Spitzer et al (1973)* which 
correspond to an excitation temperature of 'u 8o° K. In contrast 
tO'HWS we treat all of the 60 accessible Lyman band limes on an 
individual basis. Upward absorption and downward dissociation 
probabilities are taken from Allison and Dalgamo (1969) and 
Dalgarno and Stephens (1972), respectively. The latter correspond 
to a mean dissociation probability per molecule of 0.26 in an 

optically thin region, which is 'v 2/3 the older value used by HWS. 

"p" 2. 

We have considered Doppler shift parameters b = (2 v ) 2 in the 
range from 1-10 km/s, the lower value being the one used by HWS 
and characteristic of thermal broadening at 100° K. The H 2 
abundance near the surface of a cloud is very sensitive to b^, 
as may be seen from Fig. 2. Values near 4 km s _1 seem to be 

appropriate to the 0A0-3 clouds (Spitzer et al, 1973, and private 
communi cation ) . 



Our treatment of the attenuation of the UV continuum is based 
on the 0A0 selective extinction measurements, as already mentioned 
above in connection with the role of the dust grains in catalyzing 
formation. The extinction Ay required in (7) four the mean total 
cross section at X per gas atom, is obtained from the selective 
extinction measurements y/ E B y through the relation 




E 


A-V 


E 


B~V 


( 20 ) 


For the general extinction we use R - 3, in accord with the discussion 

by Bless and Savage (1972). We have adopted a constant albedo a = 1/2. 

Witt and Lillie (1973) have deduced values of a from) measurements of the 

o 

-diffuse galactic light which rise rapidly as A decreases from 2000 A 
o 

to 1500 A. It is difficult to comprehend such large values on the 

basisof reasonable grain models, and so for our present model 

calculations we have adopted a "dark grey" grain scattering picture 

with a = 1/2. The absorption cross sections a _ (X) used 

gr ahs 

here are given in Fig. 1 for the "average” extinction curve of 
Bless and Savage (1972) extended by 0A0-3 measurements (York et al, 
1973). To indicate the deviations possible, cross sections obtained 
on the same basis for £ Oph are also given. 

In line with our objective of obtaining an overall picture of 
diffuse clouds, we have used incident radiation fields characteristic 
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of some typical or average location in the interstellar medium. The 
same isotropic intensity I^( X ) characterizes the field at every 
point on the surface of the cloud, and we assume that the material 
inside the cloud has the same macroscopic velocity and Doppler 
parameter b^. These restrictions might well have to be relaxed 
in attempting an interpretation of observations toward a particular 
star. In particular, proximity of a cloud to a star ion cause a 
large increase in I Q (X), as is documented by Herbig's calculation 
(1968) of photoionization coefficients for C Oph and Jura's calcula- 
tions (1974) of fluxes for 0A0-3 stars. At the start of this work 

there were no calculations of the average galactic radiation field 

. 0 

which adequately covered the entire spectral range from 912 A to 

o 

2028 A . Thus we interpolated on Habing’s calculations (1968), 

o 

extrapolated Lillie's data (1968) below 2000 A , extrapolated 

o 

Witt and Johnson's calculation (1973) below 1330 A , and most 

o 

recently extrapolated Jura's calculation (1974) below 1125 A, 
to obtain the four curves in Fig. 2. Many of the illustrative 
calculations given below use curve 2, since Habing's field has 
been used by so many other authors. 

The intensity is attenuated in our model according to the 
formula 

« 

I(A,x) = I 0 (A) e “ T X (x) , (21) 

where the optical depth at a distance x into the cloud T.(x) 

A 

is the sum of optical depths due to absorption by grains and 



molecules. We use (21) in such a way, however, as to take into account 
approximately the fact that radiation can reach the point x from any 
direction for a finite cloud. In particular, the intensity is con- 
verted into an effective flux F(X,x) for absorption 

F( X, x) - c it I (X, x) , (22) 

with c chosen separately for the lines and for the continuum. When 
this is done, it is possible to use the method devised by the appendix 
of HWS for treating the dissociation of and to treat also the 
pressure, thermal, electrical, chemical balance relations (3) - (6). 
Values for c con ^ a nd G -^ ne were chosen on the basis of exact solu- 
tions of radiative transfer problems (without scattering) for slabs and 
spheres, the former restricted to absorption coefficients symmetric 
about the mid plane and latter to constant absorption coefficients. 
These solutions involve exponential integrals as well as simple 
exponentials, and depend on the distance from the surface measured 

by t and the size as measured by the optical depth t to the 

o 

center (X is suppressed here). For a given or a range of 

T q , one can choose c so that the difference between the absorption 

calculated from (22) and from the exact solution is minimized. For 

the Lyman band lines, for which t q » 1 as long as N > lO^cuT^, 

C line = 1 gives ver y high accuracy independent of whether a slab or 

sphere solution is used. For the continuum, where t £ 1. 

o ’ 

there are important differences between the two geometric models. 

Of course we generally do not know the geometry of diffuse clouds. 
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which might actually be quite different from either of the two simple 
cases we’ve treated. Furthermore, even if a cloud were approximately 
spherical, we would not know the line-of-sight intersection. For 

'Y- 

0A0-3 clouds with considerable H 0 abundance c .= 2 and 3 

2 cont 

respectively for slab and sphere. Using fixed values for c 
leads to errors in the radiation transfer which are typically of 
the order of 10 %, but as much as 25 % near the surface- These are 
considerably less than the differences in the solution between the 
extreme geometrical forms which might be possible for interstellar 
clouds. In most of our examples discussed below, we use c con ^ = 3.0, 
but also consider the effects of varying c con ^ from 1-4, i.e. the 

range from very thick to very thin clouds (see Fig. 3)* 

( INSERT PAGE 21a . ) 

D. Heating and Cooling 

As mentioned earlier, CII is the most important cooling agent 
for diffuse clouds. Because the theory cf the cooling by CII and the 
other contributing coolants is extremely familiar (see for example 
the review by Dalgarno and McCray, 1972), we omit any detailed dis- 
cussion of this process. We should mention, however, that the 
cooling rates are proportional to the relevant abundances so that 
depletion reduces them, thus tending to produce higher temperatures 
(Field et al, 1969; Meszaros, 1972). Further, the column densities 
of the heavy atoms and ions in diffuse clouds are small enough so 
that there is negligible trapping of cooling radiation. 

There is much less certainty as to the means by which the 
interstellar clouds are heated. We have recently published a 
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( IMSERT ) 

It is, of course, possible to treat the radiation problem for 
spherical geometry by solving numerically for I(A,x,Q) in the manner of 
HWS. Exact solutions can also be found for a slab. We have chosen to 
work with the approximate solution (21) and (22) (a corrected semi-infinite 
slab) because it allows us to calculate models for clouds of varying 
total column density in one computational run. In addition, we have 
chosen to investigate the effects of varying many of the input parameters, 
which is only feasible with the efficient numerical techniques associated 
with our sample model. 
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detailed study of the heating and ionization of an gas by cosmic ray 
protons and hy soft X-rays (Glassgold and Langer, 1973b) which can 
easily be applied to the presently considered diffuse clouds containing 
a mixture of H and First of all, the X-rays can not penetrate 

beyond a thin surface layer of the cloud since the column density 
corresponding to an absorption optical depth of unity is 


N- • « 2 x 10 19 cm** 2 ( 

1 XR 


100 eV 


E 


XR 


Within this distance, the observed soft X-ray flux down to = 100 eV 

would produce an ionisation rate g m 3 x 10" s’" and a heating 

XK 

efficiency of ^ 1/3. Although these X-rays may have some observational 

consequences, we ignore them in our model since they are inoperative 

""throughout the bulk of the clouds under consideration, which ha^ ^ column 

20 -2 

densities to the center ranging from N - (l ~ 10) x 10 cm . Despite 
qualitatively different features for the interaction of secondary 
electrons with H^, the calculated heating and ionization rates are 
not very mu^h different than for H. The heating Q and number of 
electrons produced per primary ionization event <J> for a 2-MeV 
proton and an electron concentration x - 10 ^ are: qC^) * 15 eV, 

Q(H) - 8,0 eV, and <J> (H^) a 1.3, 4> (H) - 1.5. As will be seen 

from the results given in Section III, cosmic ray heating alone can 
not produce the observed temperatures of interstellar clouds even 

with substantial coolant depletion unless unreasonably large values 

> -15 -1 

of the cosmic ray ionization C ^ 10 s are assumed. We have 
therefore been led to consider a number of other heating mechanisms 
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to see whether either singly or together they miglii not contribute 
significant heating. 

(l) Dissociative Heating 

Stephens and Dalgamo (1973) have suggested ttaat the 2 H atoms 
from the photo dissociation of may cause significant heating of 
the gas. The heating rate nay be written as 

r diss n n ^ H 2 ^ \ A E h / diss 9 


where rj is the photo-dissociation rate discussed above, and 
is the mean heating per dissociation: 


\ h 


diss 



diss 


Dp A E 

v'JJ’ OJ, v* J 1 v' 


l p 

v‘JJ r OJ, v’J 1 


(24) 


In (24), n 


OJ,v'J' 


is the rate per for dissociation via B v r J’ XOJ 


and A E^, is the corresponding mean kinetic energy released, calculated 
by Stephens and Dalgamo (1973). Eq. (23) is based on 100$ efficiency 
for conversion of H-atom kinetic energy into gas heating. For a very 


/ 

diffuse cloud / A E 

\ h . 


diss 


* 0.42 eV. By using (S) we find for 


p » R^n, i.e. the outer portions of a cloud. 


( 0 ) 

diss 


= Rl n 2 <^A 


diss 


- 6.72 x 10~^ ergs cm^s ^ r? x ( 


H. 


. n -17 3 -1 
10 cm s 


) (25) 
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r<}iss wil1 s l° w iy decrease going into the cloud. 

(2) Grain Photo electron Heating 

Watson (1972) has advocated reconsideration of grain photoelectric 
heating (Spitzer, 1948) "because the photoelectron yieM above 10 eV is 
now known to be considerably larger for some representative materials 
than the values used in earlier estimates. This point must still 
remain uncertain since the surface properties of interstellar grains 
are unknown. Watson (1974) has also noted that certain ice mantles, 
which are frequently postulated in grain models, may The inefficient 
photoelectron producers. On the other hand, irregularities, impurities, 
and small grain size, all help to increase the efficiency. In any case 
it now seems appropriate to consider the effects of grain heating on 
the thermal properties of interstellar clouds, 

/ Following Watson (1972) we adopt an extremely simple formula for 
grain photoelectron heating in view of the numerous uncertainties in 
the physical properties of the emission process, 


gr ph el 


= y 


A 1 

\ AE h> !) < * a- -_(*> To 

* err* ^ 


gr ^ o 
912A 


) , ^ A ) j C ■ 1 

gr abs cont 




( 26 ) 


The cut-off A^ is the wave-length below which the photoelectron 
yield assumes a relatively constant value according to Watson (1974). 
Both the yield and the photo electron energy have been extracted 

from integrals as the mean values y and A I \ . In particular, 

° / v X 'gr 

A^ = 1240 A, E^\ =1,8 eV, and very uncertainly 

/ gr 
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0.03 < 3 < 0.2. The quantities in the integral (26) have been 

discussed above. For a very diffuse cloud without attenuation, we 
can estimate (26) approximately in terms of a mean grain absorption 
cross section and a mean radiation field* 

r ^ - 1.25 x 10 2 ^ergs s ^ n 

gr ph el 

x (y / 0.1 )( 1/4 x 10~ 9 (cm 2 s str fr)" 1 )(d abs /10“ 2i cm 2 ) 

( 27 ) 

Equations (26) and (27) are based on 100$ heating efficiency for the 
heating of electrons with mean energy v 2 eV (Glassgold and Langer, 
1973b). Possib3.e effects of grain charge have also been ignored since 
the best estimates of the grain potential energy U (Feuerbaeher et al, 

1973) indicate that |U| « • 

(3) Formation Heating 

Spitzer and Cochran (1973) have proposed that newly formed 
molecules have considerable internal and translational energy in 
order to account for the increase in velocity dispersion with rotational 
quantum number J observed with 0A0-3. They have also suggested that 
the kinetic energy, which they estimate to be ^ 3 eV, can heat the 
gas. We have included this new heating mechanism in our program as 

r form * ^l n ^ \ y form 

= 7.16 x 10 ^ evgs erne's ^ (1 - f ) 

x ( / A E h \ /4.48 eV ) ( R 1 /10 -17 cm 3 s _1 ) . 

' / form. 


( 28 ) 
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The relative importance of the heating mechanisms investigated here 
depends sensitively on parameters which are not directly observable 

and which are also difficult to estimate theoretically- The only one 

\ 

which can be calculated with confidence is dissociative heating 
(Stephens and Dalgamo, 1973). Thus we will vary the parameters in 
each heating mechanism over a fairly large range in order to under- 
stand its potential for heating. Finally, for comparison we give the 
cosmic ray heating for a mixture of H, He, and (see Appendix B). 

r cR * (1 + ?(He) ) S (H)n(8 eV + f 7.6 eV) 

( i 

III. Results 

We present our results through a series of figures and accompanying 
comments. Most of the figures are in the form of fractional abundance 
F and neutral C column density KlfciJ vs. total hydrogen abundance N, as in 
Figs. 2-4? 6-10. F is defined as 

F * 2N [hJ/n * 2N [h 2 ]/(n[hJ + 2n[h 2 J). (30). 

In this type of figure, there are always two common curves (one for F and 
one for H M ) against which the effects of changing individual model 
parameters can be gauged. The parameters for this "standard” run are 
given in Table II. Most of the reaction rate constants are listed in 
Table A of Appendix A. The parameters in the standard run have not been 
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determined by any systematic analysis, and we do not expect them to provide 

quantitative agreement with any particular cloud. On the other hand, this 

standard run is qualitatively similar to those we might use in attempting 

to obtain detailed fits to observations for diffuse clouds. 

In Fig. 2. we see that changing the Doppler shift parameter b^ has 

a large effect at small column densities N < 1 x 10 cm - , but essentially 

no effect on other constituents such as C. We note that the large variation 

in F occurs in the range N * (10 ' - 10 ) cm , which is where the self- 

shielding of the dissociating Lyman band line radiation of molecules 

takes place. The narrower the lines, i.e. the smaller b^, the more rapidly 

(as a function of N) do the Doppler cores of the lines become optically 

thick and newly formed molecules protected from being dissociated. 

Hence the regular increase in F as a decreasing function of b^ in Fig, 2. 

In the transition region in Fig. 2, the strongest Lyman line has its 

la 

Doppler core completely absorbed for N = 10 cm , i.e. the optical 
depth is e for this N at the frequency where the Doppler and radiative 

line profiles are equal. By the same criterion, essentially all Lyman 

. 20 —2 20 -2 
band line cores are absorbed when N * 1 x 10 3 m . Beyond N > 10 cm , 

the absorption occurs in the radiative wings. In this region the variation 

in abundance can be obtained in closed form if line photo-dissociation 

is the only destruction mechanism (Solomon and WicKramasinghe, 1969; 

Jura, 1974). 

We have compared our calculations of abundances with those of 
Holleribach, Werner, and Salpeter (1971). When account is taken of numerous 
differences in methods of calculation, e.g. in radiation field, geometry, 
grain absorption, transition probabilities, etc., satisfactory agreement 
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is obtained with the results presented in their Figs. 1 and 3- As noted 
by Spitser et al (1973)? the abundances of observed with Copernicus 
are in qualitative agreement with the predictions of Hollenbach et al 
and thus with the type of curve presented in Fig. 2. 

At the surface of a diffuse cloud, (11 ) applies with a dissociation 
rate p^ >> R^n unaffected by attenuation, i.e. 


R i n o 


(31) 


For the standard run in Fig. 2 (b^ - 4 Lm/s), n Q = 32 cnf^, p Q = 0.44 x lO -1 ^ 1 , 
R^ = 3 x 10 17 s so that f^ - 2 x 10 ^ . This value is actually four times 
the value calculated in our program for N = 0 and plotted in the figures. 

This is because it is based on the value c 2_-j_ ne = 4, i.e. on the assumption 
that. .the surface of a diffuse cloud is fully illuminated by the interstellar 
radiation field. Large departures from f = 2 x 10 can. be obtained for 
different choices of the parameters as we have argued earlier (Glassgold 
and Langer, 1973c). Jura (1974) has recently discussed the effects of 
changing p^ which arise from the proximity of a cloud to a star and from 
shielding of the dissociating radiation by other clouds. 

The effect of varying is given in Fig. 3; c eont = 4 corresponds 

to an optically thin cloud whereas c = 1 would be appropriate to a 
very large optical depth. Only the extreme limits of the slight variation 
with F are given in Fig. 3. The variation of N jjCl] is much larger. Re- 
calling the discussion near the end of Sec. II. C, we note that exact 
solutions suggest - 2 for a slab and - 3 for a sphere for 

optical depths to cloud center t q * 1. For a given geometry, errors 
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made by using a constant ® con -j. are 
large as 25% near the cloud surface. 

The effect of varying the pressure p Q is illustrated in Fig. 4« Both 
F end N [Clj increase with p q because the density n increases with p Q . The 
temperature profiles for the four calculations of Fig. 4 are shown in Fig. 
5(b). For reasons to be discussed later, the temperature varies only slight- 
ly with N for the type of model defined by the parameters in Table 2. The 
mean temperature also changes little when substantial changes are made in 
p^ because T is controlled by Boltzmann factors in the cooling transitions 
in (4), i.e. only small changes in T are needed to accommodate large changes 
in other quantities in the balance equations. Thus the major effects of 
changing p Q can be traced to a change in density. The mean temperature is 
usually given for each calculation as a label located near the F vs. N curve. 
The density profiles for H and H are plotted in Fig. 5(a). Most of the 
increase in n is associated with conversion of H to at constant T. If 
we set n(H) + n(H^) * (l - — ) n s P q AT and replace T by T, then the 
variation in n can be approximately accounted for by 


typically 10% in N L CI 1’ but may be as 


n 


p At 

o 

T-f/2 


The. curves in Fig. 6 show the effect of varying the cosmic ray flux 
parameterized by the primary ionization C^(H) P er H atom. The actual 
ionization rate is increased through the action of the primary and secondary 
electrons, and depends on the composition of the medium. The rates for H, 

H and He are listed in Table A , and the effective rate for a mixture of H, 
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H 2 , and He is discussed in Appendix B. The range considered for £ goes 
from the large value ( 10 15 s 1 ) used a few years ago in the 2-phase model 
to a value ( 10 s ) about a factor of 2 smaller than the extrapolation 
by Spitzer and Tomasto (1968) f r0 m the observed flux above 1 6eV. The 
main effect of changing £p in Fig. 6 is to change the electron density, 
in accord with (17). This is the reason for the large changes in N/fciJ , 
as dictated by the C balance based on (12). The mean electron densities 
are included as an extra label on the N fcij curves in Fig. 6. Another 
aspect of Fig. 6 is that the temperature is altered because a significant 
part of the cooling comes from electron collisions with 0 + ions. This 
then has an effect on the density through the gas law (3), and consequent- 
ly on the H 2 abundance. When £p approaches 10" 1 V 1 , another factor 
enters into the determination of T, heating by the cosmic rays. The 
mean temperature for Cp = ICf 1 V 1 actually .increases relative to the 
standard run, despite the enhanced cooling arising from the factor of 
4 increase in n e> In this case cosmic ray heating becomes the largest 

of the heating mechanisms, i.e. r p = (0.5- 0.7) r 

* tot 

Another aspect of the ?p = lO^V 1 curve in Fig. 6 is the large 
reduction in H 2 abundance for N > 2 x 10 20 cif 2 , which occurs because 
cosmic ray ionization of H 2 has become the major mechanism for H 
destruction. Table III presents the total column density at which cosmic 
ray destruction is 1/2 the photo destruction rate. For the smaller 
values of £p io 16 s ' 1 , cosmic ray destruction is not very important 
for N < 10 21 cm so that the H 2 abundance ratio F can not. provide much 

more than a very rough upper limit to £p . Furthermore, F is sensitively 
affected by other as yet undetermined parameters such as the W field I 
and the formation rate ^ (see Figs. 8 and 9). Perhaps the most important 
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conclusion here is that the abundance of atomic H at large column densities 
> 21-2 

N ^ 10 cm is influenced by cosmic rays, as suggested by Solomon and 
Werner (1971). It must be realized, however, that a quantitative evaluation 
of n(H), i.e. to better than a factor of 2, is impossible unless the 
branching between the dissociative recombination channels for (3H vs. 
h 2 + H ) is known. 

. In Fig. 7 the effect of changing the depletion of C is portrayed. 


Besides the direct obvious increase of N 


[ci] wi 


with the relative abundance 


5(C), increasing 5(C) lowers the temperature by increasing the cooling 
rate. This then increases the density, which generally promotes the 
formation of and explains the upward shifting of the F curves with 
5(C). The increase in density further amplifies the abundance of CI since 
both the total C density and the electron density are increased. The 
curve for the solar abundance 5(C) -4 x 10~^ seems to be well above 
the preliminary Copernicus determinations of N W (Morton et al, 1973). 
The temperature Is also significantly lower in this case than for depletion 
factors in the range from 4-&. 

The curves in Fig. 8 exhibit the effect of changing the radiation 
field. The labels 1^ - 1^ refer to the four radiation fields given in 
Fig. 1 and discussed in Section II. The main trend of the curves can 


again be understood qualitatively in terms of the thermal balance appropriate 
to each field, when we recall that the sequence I- , I 0 , I , I correspond 
to increasing values for the total integrated UV flux. Since the main heating 
mechanisms in the model are proportional to the flux, the temperatures in- 
crease and the densities decrease as we go through this sequence. The 
calculated changes in densities are actually larger, however, because of 
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the varying degrees of conversion to molecular hydrogen. The four radiation 

fields have different shapes, and their relative abilities for dissociation 

o 

are different than for heating. (The fluxes at 1000 A are ^1.5, 2.4, 4.0, 

_ c} 2 _ 1 

and 5.0 x 10 photons/cm s Hz for the four fields 1^, I^, I y and 1^, 
respectively. ) Furthermore, the relationship between dissociating intensity 
and gas density is non-linear, so that it is difficult to sort out the causes 
of the changes in these curves. 

Some effects of varying the rate constant for formation on grains 
are shown in Fig. 9. The variation in R^ by a factor of 50 leads to the 

very large variations depicted for F, but the differences tend to he 

>20-2 I 1 

compressed once N 'u 10 cm . The spreading apart of the N I Cl I curves with 

increasing N arises in a complicated way from several factors which lead to 

small decreases in n and n(CIl) with increasing R . These changes arise 

e i 

from an increase in temperature associated with a significant reduction in 

20 

CIX cooling via excitation by H atoms. For example, n(H) at N = 5 x 10 

—3 — IS 3 — 1 _ o 

decreases from 'u 30 cm for - 5 x 10 cm's to ^ 4 cm ^ for 

—16 3 —1 

= 10 cm's . We are assuming that collisions are ineffective 
in exciting CII (Glassgold and Larger, 1973a.). 

In the above figures we have attempted to show the effects of changing 
individual model parameters. These changes can be very substantial, but 

A 

changes in one parameter can often be compensated by changes in another, and 
this will surely make it difficult to obtain a unique set of parameters. 

In Fig. 10, the pressure and the radiation field have been simultaneously 
increased by factors of 2 and of 5 in the two sets of solid curves. The 
change in radiation field is achieved by simple scaling of the Habing 
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field I^. Increases by these amounts might arise from the proximity of a cloud 
to a hot star (Herbig, 1968; Jura, 1974). The exact amount will depend on 
distance, and one side of the cloud will be illuminated more than the other — 
an effect not included in the present calculations. The curves in Fig. 10 
have been plotted on an expanded scale to illustrate how little the relative 
H abundance F is changed. This result is not too surprising in view of 
the fact, already noted in connection with Fig. that the main effect 
of varying p Q is to change the density. According to (ll), the fractional 
abundance is determined by the ratio of UV dissociation rate to density, 
and in Fig. 10 this rat5.o is nearly constant. The relatively small increases 
in temperature have a complicated explanation in the changes in the n-dependence 
of the heating and cooling rates and in the relative importance of e and H 
excitation of CXI. The electron density follows (17) and, consistent with 
(12) for n(CIl) » n(CI), leads to the decreases in Fig. 10 which are 
roughly proportional to V2~ and / 5 . 

Model calculations for large p^ and n are of interest in connection 

3 -3 

with the high densities ^ 10 cm suggested for the £ Oph cloud by Herbig 
(1968) and by Black and Dalgarno (1973b). The latter authors have suggested 
that increasing the radiation field and the density may account for the 
observed rotational state populations (Spitzer and Cochran, 1973) in 
the context of a collisional-radiative steady state. Although the same 
range of fractional abundance F can be achieved, such models imply 
substantial changes in other constituent densities. Moreover, the 
very large kinetic pressures required (10^ - 10^ °K cm require 
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reconsideration of the stability of such cloud masses. 

The temperature profiles for the standard parameter set of Table II 
were given in Fig. 5(b) along with those for different pressures. The 
standard run yields a mean temperature of T a, 63° K, which is somewhat 
lower than the ^ 80° K deduced by Spitzer and Cochran (1973) for 0A0-3 
clouds. When we reexamine the T labels in Figs. 4-10, however, we see 
that many of the parameters affect the temperature so that only relative- 
ly slight adjustments would be required to boost T to 80° K. In addition 
to varying p Q , £, 1^, R , and especially £(C), we can also vary somewhat 
the strength of the heating mechanisms. 

The temperature profiles in Fig. 5(b) are fairly typical of our 
model calculations for diffuse clouds. Deviations from T tend to be 

pi -2 

< 10% for N < 10 cm . The approximate isothermality stems from the 

near constancy of the important heating rates in Eqs. (25) - (29). The 

component heating rates are shown in Fig. 11(a) for the standard run and 

in Fig. 11(b) for the case in which the photoelectron efficiency has been 

increased to y = 0.2 (f - 80° K in this case). Dissociative heating 

accounts for only a few percent of the heating and has been omitted from 

the figures. Cosmic ray heating is minor but not negligible when 

r (H) - 10 °s , but it would become negligible if £ were reduced 
*p P 

-17 -1 

to 10 s or less. In the standard run, grain photoelectron heating 
and formation heating are comparable, accounting on the average for 
45$ and 40$ of the heating, respectively. If we had to rely on just one 
of these two heating mechanisms, then T would drop to 40° K if no 
other changes were made. But we might then also consider increasing 
the strength of the single source as well as making slight changes in 
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parameters such as p , I, R, and £(C). However, we would he hard pressed 

to achieve T = 80° K with just formation heating because the mean heating 

yA E \ associated with a newly formed E n molecule is limited by 

X V form 

observations (Spitzer and Cochran, 1973). In the standard run, we already 
used two-thirds of the maximum allowed value of 4.5 eV, and it is really 
difficult to comprehend unit efficiency in this process. Grain photoelectron 
heating seems more promising to us in that efficiencies y in the 10-20$ range 
can provide T of 00-100° K. Thus T = 80° K for the calculation in Fig. 11(b), 
where the only change from Table II is y =0.2. To illustrate the effect 
of several rather not unreasonable changes, we found that when p - 1500, 

£ = 5 x 10 17 s \ £(C) = 5 x 10 J , = 5 x 10 ^ 7 em^s \ T reached 

140° K and 105° K, at the surface and interior, respectively, for an 
average of T ~ 110° K. Of course the other effects of such changes would 
have to be confronted with relevant observations before drawing any con- 
clusion other than that a fair range of temperature is possible in this model. 

We have also applied our model to the Ca ionization equilibrium, 
making specific use of (12), ( 13 ), and Table I. Our calculations differ 
from earlier ones (e.g. Herbig, I960; Brown, 1972; Pottasch, 1972; White, 

1973) in that our photoionization rate is at least a factor of 2 larger, 
reflecting more recent cross section measurements (Carter et al, 1971; 

Mcllrath and Sandeman, 1972). In addition, our model is generally more 
complete; it includes, for example, inhomogeneous electron and hydrogen 
densities. The calculated density ratios n( Call )/n( Cal) decrease 
typically by a factor of 2 when N varies from 0 to 0 x 10 cm . This 
is somewhat less than expected solely on the basis of radiation transfer 
considerations, and reflects the effects of decreasing electron density 
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and Call! density going into the cloud. Density inhomogeneities do not 

cause appreciable differences between N gall] /N {fla j and n(c a II)/n(Cal) 

20 

until N > 5 x 10 cm, where they are 'v* 10$. 

Our calculated N ^CallJ /M [calj ratios are typically of order 

'i 

of magnitude 10, whereas the values deduced from interstellar absorption 
lines (Pottasch, 1972; White, 1973) are closer to 10 . Using a sensitive 
synthesis technique, Bortolot, Shulman, and Thaddeus (1974) have reported 
a ratio very close to 100 for X, Oph. Although the calculated ratios can 
be reduced somewhat, it is hard to obtain a. factor of 10 reduction with- 
out running into other difficulties, such as low temperatures or requiring 
high cosmic ray fluxes, weaker UV fields than Habing's, and C abundances 
near solar. 

Another way of describing this situation is that the electron densities 
deduced from measured Cal and Call interstellar absorption, lines are about a 
factor of 5 to 10 larger than given by our model. The values suggested seem 
to be "v 0.1 em~^ (e.g. White, 1973). These determinations are based on a 
single cloud component picture in which the various factors in the ionization 
balance equation (12) are replaced by mean values, and densities by column 
densities. Such a simplified picture may not be appropriate, however, and 
important contributions to the relevant column densities may come from the 
intercloud medium and/or the immediate neighborhood of the star. For 
example, the recombination coefficient for Call may be at least as big in 
a hot intercloud medium (as in a cool cloud) because of dielectronic 
recombination (Pottasch, 1972). If the electron densities are about the 
same in the intercloud and cloud phases, as is often stated (see the review 
by Dalgarno and McCray, 1972), then the intercloud region may play an 
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important role in the analysis of Call to Cal column density determinations. 

The Copernicus observations (Morton et al, 1973) only determine N^CIlJ 

to within an order of magnitude for the four stars, £ Per*, a Cam, £ Oph, 

and y Ara. Since most of the gas-phase C is in CII, these measurements 

correspond to depletion relative to the solar abundance by a factor in the 

range from 2 to 20. The measurements of N [cij are more precise, and the 

average for five stars (the above mentioned plus X OriA) seems to be about 

'u 3 x 10 ^ N. The values of N obtained from our standard run seem 

to be somewhat smaller by a factor of 'u 1| to 2. We do not believe that 

this difference should be taien too seriously, since the standard parameters 

were rather arbitrarily selected. The calculated values of N LciJ are quite 

sensitive to the choices in certain parameters as shown in the above figures. 

Larger values for N 03 could be obtained, for example, by reducing c coir j.* 

I, and R , or by increasing p , £(C), and £ . These effects can be 
1 op 

understood in terms of the basic balance equation (12), which we rewrite 
for C with the approximation n(CI )/n(CII ) - n a( CII )/G( Cl ) < < 1, as 


n( Cl ) 


n a(CII ) 


€(C) n 


Additional simplifications would usually be made in applying (32) to 

measured column densities, which could cause substantial errors. For 

example, N £*Clj /N will become substantially smaller than n(Cl)/n for 
20 -2 

N > 3 x 10 cm simply because of density inhomogeneities. With caution 
we can interpret the 0A0-3 observations to imply an average n(CI)/n 'u 
3 x 10 ^ and then use (32) to estimate n^£ /G 0 (CI) ^ 2.6 x lO^cm^s \ 
where G q (CI) is the Cl photoionization rate at the surface of a diffuse 



cloud. In other words, in this interpretation the present measured values 

of » [Cl] determine only the product of three relatively poorly known 

quantities. If we use Jura’s field (G 2.5 x lCT 1 ^" 1 ), then 

o 

n e ^ °*°3 cm ^ and ^ 0.003 cm - ' 5 for £(C) = 2 x 10“^ and 2 x 10“ 5 , 

respectively, corresponding to depletion factors of 2 and 20. Although 

this argument is very crude, it suggests that our model calculations may 

be able to explain the 0A0-3 Cl measurements with fairly small electron 

densities and cosmic ray ionization rates £ (H) < 10” 16 s -1 . This situation 

is almost the reverse of that discussed above for Ca. The Ca observations 

have been interpreted by White (1973) to imply electron densities much 

larger than are obtainable with the ’'intermediate” cosmic ray ionization 
—16 —1 

rate C (H) * 10 s . Our own. preliminary qualitative conclusion 
from the 0A0-3 Cl data is that H ) , u 10 "^s consistent with the 
crude upper limit deduced from relative abundances F in the thicker 
0AQ-3 clouds. The same conclusion has been reached by Black arid 
Dalgarno (1973a), Jura (1974), and O'Donnell and Watson (1974) by 


consideration of the HD column densities measured by Spitzer et al 
(1973). These authors obtain values for which generally lie in 
the range lCf 17 s -1 - 10“ 16 s _1 . 


IV. Conclusion 

The results presented in Section III suggest that the model defined 
in Section II may be useful in interpreting interstellar absorption lines 
arising from diffuse clouds. Although detailed analyses of individual 
clouds must still be carried out, the model does seem to be capable of 
reproducing the general trend of H 2 fractional abundances and kinetic 
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temperatures observed by 0A0-3. More quantitative tests of the model should 
be carried out as more extensive data becomes available, e.g. of the type 
recently reported by Spitzer, Cochran, and Hirshfeld (1974) Many 

applications, improvements, and extensions of this model are suggested by 
our preliminary discussion in Section III. We are currently engaged in a 
more complete investigation of the roles of other molecules which would - s . 
be applicable to both 0A0-3 clouds and to thicker ones. 

This work has been supported by Grant NGR-33-016-196 from the National 
Aeronautics and Space Administration. 
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Appendix A. Ionization Balance - . 

The reactions required for treating the electrical balance equation (5) 
are listed in Table A. Our treatment is restricted to molecular hydrogen 
clouds in which the ionization of C is still dominated by the interstellar 
UV field, i.e. N < 10 cm . Oppenheirner and Dalgamo (1974) have recently 
considered the ionization of denser interstellar clouds. Previous treatments 
of this problem under more simple conditions have been given by Werner (1970), 
Solomon and Werner (1971) and de Jong (1972). The total rate for A6 is much 
larger than was given by Solomon and Werner since they ignored electronic 
dissociation of (Glassgold and Langer, 1973b). 

The balance equations arising from the reactions in Table A are: 

H 2 : R]_n(H)n + ') + + H)n e n(H 3 + )= ( ?2 + n)n(H,,) + K 2 n(H 2 )n(H 2 + 

H 2 : "* = r ^ H 2 ' Cqr-(H) + K 2 n(H 2 ) + B 2 n e 3 

H 3 +: K 2 r ‘( H 2 M H 2 + ) = H + ) 

H + : ^(H) - ?2 (H + )n(H 2 ) + n(H 2 + )n(H) = a(H*)n e n(H + ) 

He + : £(He)n(HeI) * a(He + )n n(Hell) 

G 

C + : Ct(CI) n(CI) = a( C + )n g n( CII ) 

with 

C 2 = ? 2 (H 2 +) + ? 2 (H+) + S 2 ( 2h ) 

= ®3 (H 2 + H) + 6 3 (3H) 

In addition, the total densities for H, He, and C are: 
n = n(H) + n(H + ) + 2n(H 2 ) + 2n(H ? *) + 3n(H, + ) 

n(He) = n(Hel) + n(Hell) = £(He)n 


) (A17) 
(A18) 
(A19) 
(A20) 
( A21 ) 
( A22 ) 

(A 23) 

( A24 ) 

(A25) 

(A26) 
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n(C) = n( Cl) + n(CII) = £(C)n (A27) 

and 

n = n(H + ) + n(H + ) + n(H 0 *) + n(HeII) + n(CII) (A28) 

The effects of the charge exchange reactions (Field and Steigman, 1972) 

0 + + H t 0 + H + + 0.002 eV 

have not "been included in the above balance equations. The potential importance 

of these reactions for n(H ) and n has been emphasized by Watson (1973), 

e 

Dalgamo, Oppenheimer, and Berry (1973), and Black and Dalgarno (1973), especially 
for dense clouds. We have included these reactions in a more elaborate theory of 
ionization than presented here, and find that 0 - H charge exchange is not 
important in determining the electron density in the relatively thin, low 
density clouds under consideration in this paper. Our quantitative estimates 
invoke significant depletion of 0 (by at least a factor of 5) suggested by 0A0-3 
observations (Morton et al, 1973) and the strong temperature dependence of the 
0 + H -* 0 + H rate. The physical reason why 0 charge exchange is not 

important for diffuse clouds is that they still have sufficient atomic H 

+ # 

for the exothermic exchange 0 + H -+ 0 + H to compete favorably with - 

■f »-{■» 

induced formation of rapidly recombining molecular ions such as OH , H n 0 ; 
etc. 

In order to express all of the densities in terms of n and the various 
rate constants in a relatively simple way, we make the approximation of 
ignoring in (Al8), that is we ignore dissociative recombination (A12) 

■i* 

as a destruction mechanism for compared to charge exchange with lip (Aid) 

+ 

and the -producing ion-molecule reaction (A1'6). This approximation 
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requires low electron density, i.e.: 

K. + (K /2 - K_ )f 

x < < X = - • - ^ - 

e 1 ■ ^2 


(A29) 


With the rate estimates of Table A and T = 100° K, x = 0.1 (l + 2/3 if), i.e. 

the condition (A29) would be easily satisfied for interstellar clouds. This 

estimate is based on theoretical estimates for and $ . The latter is 
2 

^ 10 smaller than many measured dissociative recombination rates for other 
species. However, we expect our approximation to be accurate even if is 
as large as 10 . 


With this approximation, the above system (Al7)-(A28) can be solved as 

follows. First, the densities for H 9 , , and H 0 can be obtained in terms 

of n and n : 
e 


n(H 3 + )/n(H 2 + ) = (K,,/2 ^Xf/x) (A30) 

n(H 2 + ) = 1/2 C 2 (H 2 + ) f [Kj + (K 2 /2 - q >0 (A31) 

f = 1/2 (-b + /b 2 - /ae )/a (A32) 

with 

a = (2R 1 n)(K 2 /2 - i^) + (1^ + Kg/2 BjOH)/^) C 2 (H 2 + ) (A33) 

b = ? 2 + n - k 1 ? 2 (H 2 + ) + (2^X21^ - K 2 /2) (A34) 

II 

c = (2R 1 n)K L (A35) 


Eqs. (A30) and (A3l) imply that the densities of the hydrogen molecular ions 
are negligible, e.g. when C 2 (H 2 + ) ^ 3 x 10“ l6 s~ 1 , n(H 2 + ) < f 10" 6 cnf 3 
and n(H 3 ) < f 5 x 10 cm . Then (A28) becomes 

n = n(H + ) + n(Hell) + n(CII) 

0 


(A36) 
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where the individual densities are given by 


, +x 

n(H ) = - 


a( H )n 


(A37) 


„(HeII) = -£-(feJ g(He)n 
a( He ) n 


(A3S) 


n(CII) = — 1 ?)°l\ ■■ (1 + _2l£lL \ 

a(C ) n ' ) 


n 


-1 


a( C + )n 


(A39) 


and is an effective ionization rate for hydrogen, 


C H “ {h + 1/2 C 2 (H 2 J + d ^72“”^ )f + 1/2 . c 2 (H+ ) f * 

We have assumed that n(Hell) < < n(Hel). When (A37)-(A39) are put into (A36), 
a cubic- equation results; we have solved the system (A35)~(MO) by iteration. 
When n(CIl) > > n(Cl), the equation for n reduces to the quadratic 


(MO) 




(Ml) 


where r,' is defined as 


-j. 

V = 5(He) c( He) . 

“ a(He’ ) 


(M2) 


If we assume that all atoms with ionization potentials < 13.6 eV are predominantly 

doubly ionized, then £(C) in (Ml) is replaced by the total abundance of these 

HA 

heavy atoms, and (A 40 ) has the solution 


x = 1/2 
e 


u 


r + 
UiA 




4 C'/na(H ) 


-1 


(M3) 


This is a generalization of earlier results (e.g. Preston and Brown* 1970) to 


include s-^d He. 
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Appendix B. Cosmic Ray Heating of H, He Mixture s 


The heating of a gas of molecules by cosmic rays has been studied 

by Glassgold and Langer (1973b). Regions containing only atomic H had been 
studied earlier by a number of authors, e.g. Goldsmith, Habing, and Field 
(1969) ( see Glassgold and Langer, 1973b, for other references). In this 
appendix we give an approximate treatment for a mixture of H, E^, and He, 
appropriate to the central parts of diffuse 0A0-3 clouds. We also discuss 
the dissociation rate £ 2 (2H). 

We estimate the heating by following the evolution of the electrons 

produced by ionization of H, and He by cosmic rays, following closely 

our earlier work. The ionization of H and yields a primary electron Ydth a 

+ 

mean energy = 20 eV. The ion leads to a characteristic chemical heating, 

Q. (He) - 7 eV, which arises from reactions (l6.b) and (16. c). The 
ion 

associated heating rate is 


r. 

ion 


(H 0 ) ■ Q. 

2 ion 


(Hg) C p (H 2 


)n(H ? ) , 


(Bl) 


where £ (H ) is the primary ionization rate for the reaction (A4). The 
P 2 

mean primary ionization rate, £ (e~), is given by 

P 


Cle") = i f C ( H ) + (1 - f) % (H) + ?(He) t (He) , 
P P 2 p P 


(B2) 


where £( He ) Is the fractional abundance of He. The primary electrons 
coming from ionization of H, H^, and He may follow any one of these 
channels : 


e + H 2 H 2 + + e + e 2 , E^e) + E (e) = 5 eV, a = 4 x 10 17 cm 2 

-*■ + e^, E^(e) = 9 eV, a - 4 x 10 ^ 7 c.m 2 (B3) 

H + H + e. 


E (e) = 8 eV, 


„ . ---17 2 

a = 6 x 10 cm 
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e + H + H + + e + , E,(ib) + E 0 (e) = 7 eV, o * 3 x 10' _1 ' cm 2 


'1 2 


■> H* + 


(B4) 


E (e) = 9 eV, 


, 1A -17 2 

a = 4 x 10 cm 


whera E(e) gives the mean energy of an outgoing electron, and a is a 
typical mean cross-section. The mean secondary ionisation rate is given 
by 


5 S ee (e_) = 


a i (H ? )n(H 2 ) + o^ChMH) ^ ^ 

)n( h 2 ) v o;( H^T hT | c p (e } 


(B5) 


where ch.(x) i- s ionization and o^(x)the total cross-section for species 


X. 


The heating due to these secondary electrons is given hy 


z a k n(x k :)< y E:) 


r — £ ( e ) n 
sec p v „ f \ 

£ a k n( Xk ) 


(B6) 


where Q, ( E ) is the heating for electrons with energy E in channel k. When 
x > 10 4 we find the following ( approximate) form to hold 


r sec = C p ( e ") (H n(H 2 ) + 8 n(H) ) eV 


(B7) 


and that the total heating r is given by 


r cR “ (1 + £(He) ) £ (H) n(8 + 7.6 f ) e? 


(B8) 


The total dissociation rate in a mixture can he calculated in a 
similar manner, yielding 


? 2 (2H) « (0.55 - 0.1 f) C 2 (H 2 ) 


(B9) 
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We have neglected the slight dependence on composition and used a mean 
value of 0.5 C 2 (H 2 + ) in Table A. 
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Table I. 


PARAMETERS FOR THE CHEMICAL EQUILIBRIUM OF C AND Ca 


Parameter 

Units 

C.I 

Cal 

Call 

a (100°K) 

in -12 3 -1 
10 an s 

7 

7 

33 

A., 

th 

0 

A 

1101 

2028 

1044 

A th - 912 £ 

o 

A 

189 

1116 

132 

a.. 

th 

nn -l8 2 
10 cm 

11 

t 

. 3.25 

0.167 

G/a* 

- 

20 

35 

.05 


Average of measured, values (Carter et al, 1971; Mellrath and Sandeman, 1972). 


For the case of Habing's radiation field (1968), represented by the average 

O / 9 O 

flux for the 1000-2000 A band as 6.8 x 10* 4 photons/caf's A . 


Table II 


PARAMETERS FOR THE STANDARD RUN* 


Parameter 

Symbol 

Value 

Pressure 

p o 

1:000 cnf 3 k°K 

Cosmic ray ionization rate 

CJH) 

_ n -l6 -1 
10 s 

(primary for H) 

P 


Interstellar UV field 

I 

I^ - curve 2 in Fig. l( 

C abundance 

5(c) 

7 x 10 3 

Doppler shift parameter 

b D 

4 tan s' 1 

Hp formation rate constant 
(on grains) 

R 1 

Q nn -17 3 ”1 
3 x 10 cm s 

Grain total cross section 

w x) 

curve for typical star 
in Fig. 2(a) 

Radiation transfer re- 

C line 

1 

normalization factors 



°cont 

3 

Grain photo electron efficiency 

y 

0.1 

Mean formation heating 


3 eV 


* Reaction rate constants are given mainly in Table A of Appendix A. 




Sa- 


lable III 

COLUMN DENSITY AT WHICH ? (H,) = 2r)(H 2 ,I) 


C, (H) (s -1 ) N (I0 20 om -2 ) 

P 






1.9 

10 

23 


33 



53 


Table A. 

IONIZATION BALANCE REACTIONS 


Reaction Rate Reference 


Grain formation of H, 


(Al) H + H ->■ H 2 

r i 

Photodissoeiation of H 2 

(A 2) hv + H 2 + 2H ‘ 

n 

f 

Cosmic ray ionization, etc. 

(A3) p + H -* H + e + p 

h = i - 5 s 

(A4) p + H 2 + H 2 + + e + p 

C 2 (H 2 ) = 2.3 c p 

(A5) p + H 2 *^H + H + + e + p 

C 2 (H + ) = 0.02 C 2 (H. 

(A6.a)p + H 0 -> H + H + p 

£ 

(A6.b) e'+ H 2 + H + H + e 

5 2 (2H) = 0.5 C 2 (H 2 

( A7 ) p + He *+■ He + e + p 

?( He ) = 1.5 C 

Photo ionization 

(AS) hv + C + C + + e 

G(CI) 

Radiative recombination 


(A9) e + H + + H + hv 

a(H ) - * x 10 . 'T 

(Alo) e + He *> He + hv 

a(He + ) = a(H*) 

(All) e + C + ■+ C + hv 

a(C) =2 x lO^V' 

Dissociative recombination 

(A12) e + h 2 + -*r 2H 

S 2 = ( 4-8 )T -2 x 10' 

(A13) e + H 3 + -+■ H 2 + H 

S,(H 2 + H) 

(A14) e + H 3 + + 3H 

e 3 (3H) 


Sec. II. B of text 
and Table II 


Sec. II. C of text 


Sec. II. B of text 
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3 * 


= 8 3 (H 2 + H) + e 3 (3H) 
= 4 X 10“ 6 T“*cmVs 


4 



Table A» ( eontd. ) 


Reaction 


Rate 


Reference 


Charge exchange 

(A15) H 2 + + H H 2 + H* 


* 6 x Kf^cmV 1 


Ion-Molecule reaction 

(A16) H 2 + + H 2 E* + H 


K 2 « 2 x lO^cmV 1 


r (H) is the primary ionization rate for atomic H per H atom. The mean - energy of the primar; 
electrons in A3-A7 is 20 eV. The rates are slightly composition dependent, but 
usually by no more than 10-20$. 

Theoretical estimate. In view of the fact that many measured dissociative recombina- 
tion rates are much larger, e.g. a, 10 2 $ 2 , we have considered this possibility 
for B 0 in the text of Appendix A. 

References: 

1. Glassgold and Larger (1973b) 

2. Spitzer (1968) 

3. Dubrovskii and OVedkov (1967) 

4* Reu, Johnson, and Biondi (1973) 

5. deJong (1972) 

6. Reynaber and Trujillo (1968) 



FIGURE CAPTIONS 


Upper part (a): Total grain cross-section per gas atom as a function of 

wavelength based on (7) and far UV extinction measurements. For the 

typical star we use o E (V) = 5.53 x 10 22 cm 2 , i.e. A^ = N( 5/3 x 10 21 cm 

For C Oph, we use A y = 3E(B-V) = 0.96, and N = 1.2 x 10 21 cirf 2 . Lower part (b) 

Various models for the UV interstellar radiation field. The solid horizontal 

bars are Habing’s upper and lower limits for the solar vicinity based on an 

albedo of 0.9 and 0, respectively (Habing, 1968). The broken horizontal . 

line is Habing’s recommended value for a typical point in the interstellar 

-17 3 ° 

medium. The broken horizontal line (l) u = 3 x 10 ergs/cm A is an 

o 

extrapolation of Lillie’s measurements above 2000 A , and represent a 
kind of lower limit to the far UV field. The dash curve (2) is an 
interpolation on Habing’s values. Curve (3) is based on Jura’s (1974) .> 
calculation, (open circles); the dash portion is an extrapolation. 

Curve (4) is based on the work of Witt and Johnson (1973) (solid points); 
the dash portion is a rather arbitrary extrapolation to shorter wave 
lengths, and represents a kind of upper limit to the far UV field. 


Effect of varying Doppler width parameter for ” small” column 
densities. The upper left curves give the fractional column density . 
of H F = 2N £h 2*J/N and the lower right the column density of Cl vs. 
N » .N [ Hj + 2N The parameters for the dash curve - 4 km/s) 

are the common parameters listed in Table 2. 


Effect of varying the factor c ^ in Eq. (22). The upper left curves 

give the fractional abundance F of and the lower curves the column 

density N [a] of Cl vs. N. The parameters for the dash curve 

(c . - 3) listed in Table 2. 

v cont ^ 



4. Effect of varying pressure p Q . The upper left curves give the fractional 
column density of H,,, F = 2N [ hJ/N and the lower right the column 
density of Cl vs. N = N [s] + 2N . The parameters for the dash 

curve are the common parameters listed in Table 2. The labels 
attached to the F curves are mean temperatures. 

Densities n(H), n(E^), and n for the standard parameters of Table 2 (top), 
and temperature T for the 4 calculations in Fig. 3 (bottom) vs. N. 

: ’ s Effect of varying the cosmic ray ionisation rate. The explanation of 
the format of these curves is the same as in Fig. 3. r (H) is the 

'p 1 

primary rate for atomic H; the total rate for H and the rates for H 

are listed in Table A of the Appendix. Note that the curve for F 
-IS -I 

when C - 10 s has not been drawn since it is practically in- 
distinguishable f?oa that for c = 10" 17 s _1 . The labels on the upper 

and lower curves are the mean temperatures and electron densities, 
respectively. 

Effect of varying the relative C abundance £(c). The explanation 
of the format of these curves is the same as in Fig. 3. The labels 
attached to the N £b I J curves are mean temperatures. The curve 
for £(C) = 4 x 10"^ corresponds to the solar abundance. 

Results for the four interstellar UV radiation fields plotted in Fig. 

1(b). The explanation of the format of these curves is the same as 
in Fig. 3. The labels attached to the F curves are mean temperatures. 

Effect of variation of the rate constant for H 2 formation on grains. 

The explanation of the format of these curves is the same as in Fig. 4. 

The N [cij curve for 1^51 lO'^cmV 1 has not been plotted because 
it is practically indistinguishable from that for R = lo" 17 Gm' 3 s“ 1 



Effect of simultaneously increasing p Q and X by a factor of 2 and of 5. 
The shape of the radiation field is that of the Habing flux I 9 in Fig. 1. 
The dash curve is for the standard parameters of Table II. The labels 
attached to the F curves at the top are the mean temperature T and the 
mean density n , and to the N curves at the bottom the mean electron 

density. Note the expanded scale relative to the previous figures. 

The total (solid line) and 3 major (dash curves) heating rates as a 
function N. The upper part (a) refers to the standard parameters of 
Table II; the lower part involves changing the photoelectron heating 
efficiency from y = 0.1 to y = 0,2. The curves for dissociative 
heating are not included because they represent only a few percent 
of the heating. 
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